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Transmission and Reflection of Electromagnetic 
Waves in the Presence of Stratified Media 

James R. Wait 

A general analysis is presented for the electromagnetic response of a plane stratified 
medium consisting of any number of parallel homogeneous layers. The solution is first 
developed for plane-wave incidence and then generalized to both cylindrical and spherical- 
wave incidence. Numerical results for interesting special cases are presented and discussed. 
The application of the results to surface-wave propagation over a stratified ground is con- 
sidered in some detail. 

1. Introduction 

In the analysis of antenna radiation and kindred problems, it is usually necessary to 
consider the inhomogeneous nature of the surrounding media. A good example is when an 
antenna is located over the surface of the earth. In many cases such as this, it is permissible 
to regard the ground as a semi-infinite medium with a plane surface. Furthermore, this half- 
space can often be idealized as a homogeneous medium. As long ago as 1907, Zenneck [l] 1 
showed that a solution of MaxwelPs equation could be obtained for such a region by matching 
certain wave functions at the interface separating the two half -spaces. His solution was 
characterized as a wave which propagated along the interface without change of waveform 
with a velocity greater than that of light and was exponentially attenuated with height above 
the interface (towards the air). 

In 1909 Sommerfeld [2] presented an elegant solution for the fields of dipole sources 
(electric and magnetic) which were located in the air over a homogeneous flat earth. In the 
final reduction of the integrals to a form suitable for computation, a small but significant error 
was made. This error was first pointed out explicitly by Norton [3], but not before many 
had accepted the earlier solution. To compound the confusion, it turned out that Zenneck's 
surface wave was very similar to Sommerf eld's solution in his 1909 paper. Since that time 
a large number of papers have appeared treating this subject. The confusion that remains 
is largely due to the differences in terminology. 

The discussion given here relates to three related aspects of this problem. The first is a 
general analysis of reflection of plane waves from a parallel stratified medium consisting of 
M homogeneous slabs, the second is the extension to a line source over the stratified medium, 
and the third is the generalization to dipoles or current elements over the stratified medium. 
A special case of the latter is when the number M becomes 1, corresponding to the Sommerfeld 
problem for a half -space. 

2. Plane- Wave Incidence 

A plane wave with a time factor exp (i<af) is incident at an angle 6 on a stratified medium 
composed of M homogeneous layers. The electric vector is in the plane of incidence (xz plane) . 
The situation is illustrated in figure 1 where the y axis is out of the paper. The electrical 
constants of the layers are <r m , e m , and ju m where the subscript m indicates the mth layer below 
the surface. 

From symmetry it can be seen that the magnetic field has only a y component and for 
the mth layer, it is a solution of the equation 

(V 2 -y 2 m )H my =0, (1) 

where 

y 2 m =ia m n m co—e m fx m o) 2 . 



1 Figures in brackets indicate the literature references at the end of this paper. 
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Fictre 1. A stratified medium consist- 
ing of M homogeneous layers. 
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The general solution is of the form 

H my =[a m e- u >»*+b m e u >n z ]e-* x , 
where 

and X can take any value. The incident field H lnc - can be written 



(2) 



H%f=H ( >e- yQOO * 9 ' z Xe- y »° ia9x . 

Therefore, in eq (2) a e~ u0z e~ iXx can be identified with H'™' if a =ff and i\=y sin 0. Con- 
sequently, b e UQZ e~ iXx is a reflected wave, and the angle of reflection is 6. 

The boundary conditions at the interface z=0, z—z x% . . 
fields should be continuous. Now since 

Emx= — (<r m +io>e m ) _1 ~^rr' 
this means that the boundary conditions can be written 



z=z m -i are that the tangential 

(3) 



/ I • \ — 1 Ott m _ \y , . . \ —I OH- m y 



(4) 



where m = 0, 1, 2, . . . M-2, M-l. 

Imposing the condition that only outgoing waves are permissible in the lowest layer 
(which is semi-infinite) it follows that b M =0. The boundary conditions then lead to 2 (M-l) 
equations that are linear in a m and b m to solve for 2 (M-l) unknowns in terms of the known 
coefficient a . The solution 2 is 



b 0= K -Z K 
a K J rZ l 



(5) 



fhere 



Z X =K, 
Z 2 =K 2 



Z 2 J rK l tanh u Y hi 
K x -\-Z 2 tanh u x hi 

Z s +K 2 tanh u 2 h 2 
K 2 +Z z tanh u 2 h 2 



(6) 



? The quantity X should not be confused with the wavelength. 
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Figure 2. Transmission line analogy for the stratified medium of M layers. 

y jr /s m +\-\- rL m tanh u m fi m 

m ~ m K m +Z m+1 tanh u m h m 



/jM-\ — %Lm-\ 



K m -\-K m _i tanh u M _ih M -i 
Km-i+K m tanh u M _ 1 h M -i 



where 



Km=-TT—> and^=(X 2 +\^. 



(Note that h M is the thickness of the mth slab). 

The quantity b /a is the ratio of the amplitude of the reflected wave to the amplitude 
of the incident wave. It is denoted by ifj to indicate that the electric field of the incident 
wave (and also that of the reflected wave) is in the plane of incidence. 

The present problem has a well-defined analogy in transmission line theory. 

In this analogy each section of the multiple tandemed line is to correspond to a slab. 
The voltage across the line is E mx and the current is H my for the mth section. The propaga- 
tion constant is u m and the characteristic or surge impedance is K m of the mth slab. The 
incident wave that comes from the left in figure 2 is given by a e~ w ° 2 , and the wave reflected 
at the junction, 2=0, is b e u °*. The input impedance of the line that is the ratio of the voltage 
to the current at 2=0, is Z x . Furthermore, the impedance at the junction z=z m is Z m+1 . With 
this analogy and a knowledge of the behavior of one-dimensional transmission lines, one could 
write down the solution of the two-dimensional reflection problem. 

Some features of the wave problem will now be discussed. The quantity Z x will play an 
important role in the following. Since 

Z\ = ILQxlHoy] Z = 0=lL Ix/Hly] 2 = 0, 

it is called the surface impedance, being the ratio of the tangential fields at the air-ground 
interface. In the case of normal incidence, 0=0 and \=0, so that u m =y m and K m =rj m where 

and 

r] m =[iiJL m co/a m + iernO)] 1 ' 2 . 

For example, in the case of a homogeneous ground (hi — >°°), 

Zi = rji, K = rj 0y 
and the reflection coefficient becomes simply 

#|| = — = j (/) 

do Vo + Vi 
207 



Another special case of considerable interest is when 6 approaches 90° corresponding to 
glancing incidence, then 

u m ^ (yi~ylY A , K m ^ (y 2 m ~ylY A /((7 m +ico€ m ) , 

which yields for the homogeneous ground 

%=&=*<!— &>&» (8) 

This special value of Z u which relates the tangential fields in the limiting case of glancing 
incidence, turns out to be an important quantity in further work. For this reason, it is de- 
noted by Z°, where the superscript v indicated that the electric field in the air is nearly ver- 
tical for |ti|>>|to|. This fact can be shown by evaluating the wave tilt, which is defined by 



TF=fkl (9) 



and is the complex ratio of the horizontal to the vertical electric field in the air just above 
the ground. It readily follows that for the general case [4] 



For a homogeneous ground 



E 0z /H 0y J z =o " To sin 6 



. v (l-4sin 2 0Y 2 

Y o sin0 rjo sinfl 



which for grazing incidence becomes 

For radio frequencies (co^lO 6 ) and moderately conducting ground (c^10~ 2 ), |y /Yi| is of the 
order of 10~ 2 , and thence PPis small and hence \E 0z \ is much greater than \E 0x \. 

To indicate in a simple way as is possible, the influence of stratification on the reflection 
of waves from the ground surface, a two-layer case will be considered. This is effected by 
letting A 2 -^°°- Furthermore, it will be assumed that |ti/to| an d ^/Yol^^ 1. Then for any 
angle of incidence 

u m =y m (l — ~^$m 2 oy c*y m (ro=l,2), 

and 

This leads to the simple relation 

Zi»QK u (13) 

where 

, (71/72) + tanh yjii 
1+ (71/72) tanh 71A 



n Win 2) -r taim y x m , 



and where Q is the correction to the characteristic impedance K x of the upper layer to account 
for the presence of the lower layer. Note then, if |7iAi|>> 1, (?gl. It can be said that the 
lower layer is not detectable when \Q\ is within 5 percent of unity. Such a condition is met 
when ( a 1 ju w ) 1/2 h^> 3 . 

It should also be noted that Q relates the wave tilts for a stratified (two-layer) ground and 
that of a homogeneous ground by W^ W Q, where 
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Wt&W]*,-.. 

An example is here quoted to illustrate the order of magnitude of the quantities involved: 

Frequency, /=co/2tt= 125 kc 
Upper layer conductivity, o-^KT 2 mho/m, 
Dielectric constant of air, e =8.854X10- 12 f/m, 
Dielectric constant of ground, e 1 = 10-e , 
Magnetic permeability, /i=47rX10~ 7 h/m. 

For these values, W =0.082/ 41.1°, for a homogeneous ground. In the case of a two-layer 
ground where hi is finite, 

W^WoQ 

= 0.082\Q\/ 4tl.l° + q, 

where q is the argument of Q expressed in degrees. For frequencies of this order, e l o)/ ( r l and 
6 2 w/<t 2 are small (in the above example ^g x = 0.0069). Therefore, n^(^iA*«)^ and y 2 ^ 
(i^fio))^. A formula suitable for computation of Q is then given hy 



Q- 



(o-j/^^+tanhV^F 



= tanh[V? F+tanh-^i/tra)**]. 
If o- 2 >>(T 1 corresponding to a highly conducting substratum, 

or if tr 2 <<o- 1 corresponding to an insulating substratum, 

Q~cot\h/iV. 



(14) 



(15) 



(16) 



In the above example, the parameter Fcan be replaced by hJ10, where h t is the thickness 
of the upper stratum in meters. The function Q and its argument q are plotted in figure 3a 
and 3b as a function of V for various values of the ratio a 2 /<ri. 
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Figure 3a. Amplitude of correction factor for a two-layer 
ground. 3b. Phase of correction factor for a two-layer 
ground. 
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3. Extension to Perpendicular Incidence 

In the preceding problem, the incident-plane wave has the electric vector contained in the 
plane of incidence (and the magnetic vector parallel to the interfaces). For this reason, it 
is termed "parallel incidence" . The other important case is when the electric vector is per- 
pendicular to the plane of incidence. This is termed "perpendicular incidence". m 

Again choosing the plane of incidence to be the (xz) plane, the incident wave now has only 
a y component of the electric field. By analogy to eq 2, the general solution is of the form 

E my =[a m e~ u m z +l m e^] e -iu 7 (17) 

where u 2 m =X 2j ry 2 m , a is the amplitude of the incident wave, b is the amplitude of the re- 
flected wave, and b M =0. In this case, 

in m <aH mx = ™ v > (18) 

so that the boundary conditions now become 



E , =E 



^m. y 






(19) 



These are transformable to the boundary conditions for parallel incidence by making the 
substitutions: E my for H my , ifi m <a for <T m -\-ie m o). Then using the previous results, the solution 
for perpendicular incidence can be written down: 

where 



for m=l, 2, 3, . . . , M— 1 and 
In the preceding 

where, as before, u m =(\ 2 -\-yffi*- 



y -\t ^m-fit^w tanh U m fl m (9\\ 

- iym N m +Y m+l tsmhu m h m ' {Zi) 

v m =n m . 



N m =^, (22) 



The quantity bo/do, which is the ratio of the amplitude of the reflected wave to the incident 
wave, is denoted by B ± . There is a similar transmission line analogy for this problem, which 
need not be pointed out. 

In analogy to the surface-impedance function Z 1} the quantity Y L is a surface admittance 
and is given by 

^i = — Ho x /Eq v ] Z=0 ——Hi X /Ei y ] z= =Q. (2 3j 

In the case of a homogeneous ground at glancing incidence (0-»x/2) , it follows that 

Y 1 =N 1 =(lM(l-yVyl) li , (24) 

which is denoted Y h , where the superscript h indicates that the electric field is horizontal in 
contrast to the near vertical electric field associated with Z v . It is interesting to note that 

r»Z*=(l- 7 ?/7?), (25) 

andfop |t5|»|t8I, F»Z*al. 
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4. Impedance Matching and Natural Oscillations in Stratified Media 

In this section, some remarks will be made concerning the nature of extreme conditions 
where the reflection coefficient on a plane stratified media becomes zero or infinite. The dis- 
cussion will be confined primarily to parallel incidence, although the results are easily carried 
over to perpendicular incidence. 

As indicated, the reflection coefficient B\\ is given by 

B|l= sSi' (26) 

where 

and Zi is the normal impedance at the interface z=0. In the case of a homogeneous half 
space (Ai— >°°), 

Z 1 =^ 1 =W(<ri+^€ 1 ) = (X 2 +7i)^/(cr 1 +^). (27) 

The condition for matching is that ffijj=0 corresponding to the absence of a reflected wave. 
This requires that 

K =Z l . (28) 

For a homogeneous half-space, the condition is simply 

K a =K u (29) 

which, when solved for X, yields 

LMo/7o-M?MJ 
Iii the case of no magnetic permeability contrast (mi = mo), the above simplifies ' A to 

and if both media are perfect dielectrics y l = ik l and 70 = ^0, 

But, since \=& sin 0, the preceding equations for a condition of matching can be employed to 
determine the angle of incidence when there shall be no reflection. In the case of the two 
dielectric half-spaces, the condition is 

si,1 * =± MW*' (33) 

or tan #=±(&i/&o). This latter equation is well known from optics, and there the angle d is 
known as the Brewster angle, and the ratio kjko is called the relative refractive index. 

It is indeed very interesting to observe that if k x is not real, (or y t not imaginary), the 
quantity X or required for matching is complex. That is, in the language of optics, the 
Brewster angle for an absorptive homogeneous medium is complex. The existence of plane 
waves with a complex angle of incidence will not be discussed at this stage, but it is interesting 
to note that the so-called Zenneck surface wave is very similar to the disturbance propagating 
along the interface, when there is a condition of wave matching. His solution is briefly 
described below. 



; In this case the \\;i\v til! I \'= 70/71, as can be deduced from eq. 11. 
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Zenneck [1] postulated that a surface wave could propagate along the interface z = 
between two homogeneous media. He assumed that in the upper medium (2<0) 

H 0y =b e +u o z e- iXz , (34) 

and in the lower medium 

H ly =a l e- u i z e-* x , (35) 

for propagation in the ^-direction. The signs prefixing u and u x were chosen to insure that the 
field decayed to zero for |2|->». The corresponding electric fields are then obtained from 
Maxwell's equations, 

E 0y =-K b e u o z e-* x , (36) 

and 

E ly =K l a l e- u ^ z e' rK \ (37) 

The boundary conditions then require that b =a x , and 

K +K x =0 or ^+ U ! . =0. 

The latter condition cannot be met for any value of X if the imaginary parts of u and u x are to 
be greater than zero. It now appears that Zenneck solved the matching equation K =K X in 
place of K = — K x . 

Returning to the general case, it is seen that 

has a pole (i. e., becomes infinite) when 

K +Z x = 0. (39) 

This is termed the "resonance condition.' ' For the homogeneous half -space, this resonance 
condition is simply K -\-K x =0, which is identical to Zenneck's equation for a surface wave. 

Although it is not possible to set up a surface wave in the sense envisaged by Zenneck, it 
is possible to retain his solutions if the disturbance in the upper medium is a wave containing 
a factor e~ u ^ in place of e +u <> z , and the boundary condition becomes K =K h which corresponds 
to the "matching" condition. 

In the case of a stratified media, it may well happen that the "resonance'' condition 
K Q J rZi=0 is satisfied. Rewriting this as 

Z x =i^, (40) 

€ CO 

and since for a surface wave u Q is to be mainly real (i. e., small imaginary component), then Z\ 
is to be mainly imaginary (in the positive sense). Such a surface characteristic is described as 
inductive. The quantity X, which describes the transverse propagation, is mainly real and 
somewhat greater than k [i.e., \ 2 =ul J rk 2 \. 

A simple example of an inductive surface is a highly conducting plane, coated with a thin 
uniform film of dielectric. From eq 6 with Z 2 =0 (or y 2 = ») 

Z X ^K X tanh^x, (41) 

where h x is the thickness of the film. For thin films, and /io=Mi> 

Z x ^K x u x h x =ifiouh x 1— p * (42) 

This can be further approximated to 

Z x =iwh x [l-(k /k x y}, (43) 
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since X is near k . Such a surface is therefore almost purely inductive. In this case, 



and 



u o =e Q ix o) 2 h=klh, 



X=(ug+*g)^=fc [l + (*oft) a ]»^fco[l + (fcoA) 2 /2], 



where h=hi(l—kl/k 2 1 ). The wave has a phase velocity in the transverse (x) direction, which 
is l—(k h) 2 /2 times that of free space, and it decays rapidly in the ^-direction. This approxi- 
mate solution is valid if k h (or k hi) is small compared to unity. 

In general, the determination of the poles in the reflection coefficient R {1 leads to transcen- 
dental equations that must be solved by numerical or graphical means. A discussion of sur- 
face waves on dielectric coated conductors without the restriction that hi be small has been 
given by Atwood [5a]. The subject has also been discussed by Goubau [5b]. 

Further remarks regarding surface waves will be made in a later section following the 
analysis of the excitation problem. 



5. Line-Source Excitation 

In problems of reflection from stratified media, it is more meaningful from a physical 
standpoint to consider the source a localized distribution of current. No physical source in 
existence produces pure plane waves, although in certain limiting cases the fields have charac- 
teristics closely akin to plane waves. An example is a line source of electric current. It was 
shown that at large distances from the source, the electric and magnetic field were mutually 
orthogonal to one another and to the direction of propagation [7]. It might then be conjec- 
tured that the reflected wave for a line source located high over a plane interface could be 
computed from the plane-wave reflection coefficient. Such is the case if certain limitations 
are imposed. 

It is the purpose of this section to extend the plane-wave solution to the case of an (infinite) 
line source of (constant) electric current located at height h above the stratified medium con- 
sidered in the previous section. The coordinate of the line source is taken as (0,y, — h), as 
indicated in figure 4. 

In view of the symmetry, the electric field only has ay component E m for any of the layers 
(the subscript y on E m is dropped for convenience). The direct, or primary, field E\ in the 
region z<CJ) (above the interface) is expressible as 



-w 



A'„[yo[* 2 +(z+« 2 ]«], 



(44) 



where K Q is a modified Bessel function. To convert the expression for E\ to a more useful 




Figure 4. Line source over a layered 
half -space. 



213 



form, it is written as a Fourier integral as follows: 

J 00 
A(\)e- m d\. (45) 

- 00 

By making use of the inverse property of the Fourier integral, it follows that 

A(X)=4- I " E p e ixx dx (46) 

ZK J -oo 

= = ?f^ r™K n [ y(> [x*+(z+hWW x dx. (47) 

This integral is evaluated to give 4 

ZttUq 

where u =(\ 2 -\-yl) V2 and where the + sign is to be employed for (z-\-h)<^0 and the — sign 
for (z+h)>0. 

The resultant field in the space 0]>2^>— h is then written in the form 



E = l ^ T r u^[e- u ^ z + h) + B L (\)e u ^ x - h) ] exp (-i\x)d\, 

47T J _„ 



(48) 



where the term containing Rjl(X) accounts for the presence of the stratified medium at 2—0 and 
is as yet unknown. Eq (48) satisfies the equation (V 2 — yl)E =0, and behaves in the proper 
manner as z->—h and #— >0. The corresponding expression for E in the region — A>2> — °° 
is identical to eq (48), except that the (— ) sign on the first exponential should be changed to a 
+ sign. E Q then decays properly to zero as z->— <». 

The integral representation for E Q has a clear physical meaning. When the symbol 
X is identified with k sin 6 or —iyo sin 0, the field Eq is a spectrum of plane waves of (complex) 
angle of incidence. The complete field E contains a spectrum of both incident and reflected 
waves [7]. The structure of the integrand in the integral representation for E is identical to 
the corresponding factor in the plane-wave solution in the previous section. In fact, the 
correspondence carries over to each of the sublayers. For example, in the present problem, 
for m=0 to M, 



Em= S- [ s -W 6 " ttw2 +^W ettm2 V^ x ' 



(49) 
and 

lH m G)£l mx =-z—t —lfj, m 0)ll mz =-^—' 

The solution of the present problem 5 is thus immediate; 

^ (X) =^=iV^X)+F^X)' (o0) 

v ^v Frnfi W+-A^(X)tanh u m h m 

ImW N m (\)+Y m ^(\)Unhu m h m ' [ol> 



where 



with m=l,2 . . ., M— 1, and Y M (\)=N M (\). 

Equation (48), along with (50) and (51), constitutes the formal solution of the problem. 
The remaining task is the evaluation of the integral in (48). Except in certain limiting cases, 
certain approximations must be made in order to obtain a useful result. 



4 Campbell and Foster [6], pair No. 917. 

8 It should be noted that eq (50) has the same form as eq (21) 
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The simplest limiting case is when the interface becomes highly conducting. Then, since 
^oo, R ± (\)** — l, and 



Eo- 



— ?/x w 



47T 



-TV' 

%J — oo 



\ e -u (z+h) e Uo( - Z ~ h) ]e~ iXX (i\. 



(52) 



Both these integrals are of the same type and may be readily evaluated in view of eq (46) 
and (47). Consequently, for £<C0, 



E =— ^/ifofrofr^ 



(53) 



It can be readily verified that (V 2 — yl)E o =0 for — £>0 and, moreover, E =0 for 2=0, which 
is the required boundary condition for a perfect conductor. It should be mentioned that an 
equally simple result is obtained for cases where mi = °° , corresponding to a rather hypothetical 
situation where the interface at z=0 is behaving as a perfect magnetic conductor. Then 
it ± (X)^ + l and the solution for E is identical to eq (53), with the second term inside the 
braces changing sign. 

When |7o^|>>1 and A>> — z, eq (53) for E can be written in the form 6 



*»= 2tt \2) ( 7o «o)» l J ' 



(54) 



where u =((\) 2 -\-yl) h , X= — iy sin 0, and where Q is the geometrical reflection point such 
that the ray IQ and the ray QP make equal angles with the normal at Q. The quantity i? 
is the length of the ray from the line source to the point of reflection Q. In fact, Q is taken as 
the new origin of the cartesian system r, ?/, z, as indicated in figure 5. 




p ( x y z) Figure 5. Geometrical interpretation for the secondary 

field. 



Equation (54) can be simplified even further b} r noting that the quantity prefixing the square 
bracket is simply the amplitude of the incident wave. When this is normalized to unity, 

Eo= [ e -u (z)_ e +u Q (z)^ e -i& i ( 55 ) 

which is readily identified as an incident and reflected plane wave with an angle of incidence 
and reflection of 0. To generalize, eq (55) for a stratified medium, the reflected wave would 
be simply modified by a reflection coefficient for an angle of incidence 0. The modified form 
would be 

E =[e-~» z + R±(l)e + ~<> z ] e~ ilx . (56) 

A suggested modified version of eq (53) to account for the stratified nature of the medium 
would then be 



(57) 



when R±(d) is the reflection coefficient for a plane-wave incident at an angle on the interface 
2—0. The angle is determined by simple geometrical considerations, as indicated above. 7 

« The first term of the asymptotic expansion of the modified Bcssel function Ko(a) is only retained, i. e., -Ko(a)^e- a (x/2a) H . 
• Note then tan "d=x/(h — z). 
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This result, which has been derived by an elementary approach, considers only the specular 
component. For smaller values of \z\ and h the situation becomes more complicated because 
surface waves within the layers may be excited. This question will not be pursued further 
at the moment. 

6. Line Source on a Homogeneous Medium 

When the line source is situated at the interface (h=0) between two media of semiinfinite 
extent with electrical properties e , ju and <r u ci, mo, respectively, the electric field E is given by 

* Jo ^1+^0 

for z>0. This is a special case of eq (48) with h=0 and z 1 =h 1 = °o . Noting that 

1 _ u l — U q_ U x — u 
u 1 +U ~u 2 1 —u 2 ~y 2 1 — yl 
it follows that, for z=0 

E Q = . 9 Mo 9X X Ui cos \xd\~ u cos \xd\ >• (59) 

nTi— To) \,Jo Jo J 

r™ e~~ u °\ z \ 

Lim I cos Xx — d\=K (y x), (60) 

z->o Jo ^o 

which is equivalent to eq (45), it is seen that 8 

(cos \x) u dX=f 7 g-— j K o (7 x)y o x= -K 1 (7 x)y x~ 1 ' (61) 

Making use of this result, the expression for the electric field becomes 



Since 



— ifJL uI 

~T(yl-yl)x' 

This formula is exact. 

At large values of the argument 



g o= .,., 2 _°lw 2 [To^i (WD -yixKi (Ti*)]. (62) 



K '^ x)0 ^f^' (63) 

and if the real part of y x x is >>1, 



= 7r(7i— 7oK 



£ ^/r?/ iM. (64) 



Furthermore, if y =^o where &o# is a real number >>1, it follows that 

T(yl-yt)x\2y x/ 

=ir(-Y{+tt)\2k x 3 ) {b0) 

It is convenient to express this in terms of the primary field, 

^^(^yv-, tee) 



2ir 



8 In eq (61) use is made of dKilda^—Ki/a—Ko, where a is the argument 70Z. 
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in the following manner: 



E _ 2 7o 



2ik 



El~{y\-y^x (Jl+kl)x 



(67) 



The vertical electric field in the interface, for the current line source in the interface, is 
given by 

(68) 



TT _TT _ 1 ^0 

H-Oz—H-u——: "3T— 



ifi o) dx 



1 



and, if the real part of y 1 x>>1, 



[2J0XK, (y x) +y 2 x*K (y x)-2y 1 xK 1 fax) -7^ 2 #o(7i*)] 



Ha- 



To/ 



0z= Hyl-yl)x 2 
which simplifies further, for the case |7o£|>>1, to 



[2K, (y Q x) +y xK (y Q x)], 



where 



H "- a i^feife)" v " aE »"»' 



7j =W>/7o==/W&o== W«o) 1/2 - 



(69) 



(70) 



(71) 



The field at very low frequencies near the line source is also of interest. For example if 

M« i, 



H BI 



I 



Jpyjr 



^2-2y 1 xK 1 (y 1 x)-(y l x) i K (y 1 x)], 



(72) 



which approaches —I/2wx as y t x approaches zero. 

The magnetic field H 0z expressed as a ratio to the primary field Hg z (= — I/2irx) is plotted 
in figure 6 as a function of |7 X a;| under the assumption that eico/cr^O. 
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Figure 6. Response of a line source on 
a homogeneous half-space. 
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7. Line Source Over a Thin Layer 

The general solution for M layer can also be specialized conveniently to the case of a line 
source at a height h over a thin sheet of thickness hi and conductivity a 1 . The following sub- 
stitutions in eq (48) are made: // w =mo, M=2, 72=70 where it is supposed that the medium 
above and below the sheet is homogeneous. The electric field for z <^0 is given by 



Eo=- 






< if (Yon 



)-K (y o r 2 ) + 



I 



> 2^o(2-») C os\xd\\ 



u e +u 



r 



(73) 



217 



where 

Uq+Ux tanh uji\ 
u e=u 1 -^ — -, — r — -r> 
Ux+Uq tanh^i/i! 

rf= (z+h) 2 +x 2 and r 2 2 = (^-/i) 2 +x 2 . (74) 

For small values of the sheet thickness hx and large values of the conductivity 0-, tanh 
Uxhx^Uxhx and eiw/(r 1 <<l, 
and therefore 

u ? l hx=(\ 2 +y 2 1 )hx^2i( L (75) 

where q=aiU o)h h l2. With these approximations, ^ e ^w +2ig and then 

Using the result appropriate for large q. 

* =iS^T(-l) w , (77) 

/ ft \m+l fi u (z-h) 

< u °«- h) =(£d —' (78) 



and noting 



it follows that, asymptotically, 9 



where use has been made of the relation 



relation 

f OT e u (z-h) 

cos \xd\ = K (yor 2 ) > (80) 

Jo ^0 

valid for (z-A)<0. 

The field E can thus be developed into an (asymptotic) series of inverse powers of qr 2 where 
the coefficients involve derivatives with respect of the modified Bess3l functions. Using the 
relations K r Q =—K x and K^—K^Kq+Kx/x, the nth derivative of K Q can be expressed in 
terms of K and K x . The first few terms of the resulting series are given below: 

Eq== -^ [K (y rx)-K (y r 2 ) + T(l r 2 ,g)], (81) 

where 

= ^%^-h)K l (7 r 2 ) + ~^yl(z-hrKo(lor,)-%r 2 (l-2 fc^tf^)} (82) 

+ terms containing l/(gr 2 ) 3 , l/(gr 2 ) 4 , etc., (83) 

where qi\~ o-i/x co/2,ir 2 /2 and z<^0. 

It can be readily verified that the field below the thin conducting sheet (z^>h) is given by 



S2=-^[r(7„r 1 , 2 )]for S >0- 



(84) 



As a check, it can be seen that E —E 2 at z= 0. Expressions for the magnetic field components 
can be found from Maxwell's equations: 

H mx =J-^ and H m ,=-J-^ for m=0,2- 

ljJL O) OZ tflQU OX 



A more rigorous derivation of this series expansion is given in reference [7]. 
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The above expansions developed in inverse powers of (qr 2 ) are particularly suitable for highly 

conductive sheets when the source or the observer are not near the sheet. Alternative ex- 
pansions can be developed for smaller values of #r 2 . These are developed as follows. At low 
frequencies where 70^ and 7^2 are small compared to unity, 7 can be set equal to zero, and it 
readily follows that 



where 



#o=-4^[log fo/rO+Totf)} (85) 

TM= r "T? QSX ^ x and p=h-z. (86) 

Jo A-M<£ 



Inserting the power-series expansion for cos \x, E^ becomes 

(—l) m x 2m r m \ 2m e-^d\ 



Since 



where 






(87) 
(88) 



it follows that 



= () (2m!) Jo \+iq 

Jo X-M? 

Ei(-i0q) = -l e -^dt< (89) 

™= ~ S L *r K£) [•* iK <-^} (90) 



where g=z(3q= (/*,— z)hi<riHQ<a/2. This is an ascending-series expansion in </. When x = 0, it 
reduces to 

T (t3) = -e i e«Ei(-ipq). (91) 

Using the latter result, it follows that the field Eq for 2 <0 and E 2 for 2>0 in the plane x — is 
given by 

^0= =g^[log J|=*[-«««-»lt^(_i| g _fc|j)], (92) 

for ,i<0, and 

£ 2 =±^? e i[z+hl «Ei(-i\z + h\q), (93) 

for s>0. 

The exponential integral Ei(-ig) with imaginary argument can be expressed in terms of 
the cosine and sine integrals as follows: 

Si(-v).=F=C«(flr)+i[|-iSi(^)], (94) 

where 

Ci(g) = -J m ^dt and Si(g)=j'^dt, (95) 

are tabulated [8] for real values of g. 

It might be mentioned that T Q (p) for #>0 can also be expressed in terms of exponential 
integrals, but then the argument is complex. For example, 

TJB) = 



1 f°V-X(0-"» 1 r m p -Mfi+ix) 



=— K e^ ff [e- <<pr ^i(— i(0— ix) g) +e ifla; £'i(— i03+tJc)g)]- (96) 

Unfortunately, the exponential integral of general complex argument is not well tabulated. 
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8. Radiation Field of the Line Source for Any Number of Layers 

The evaluation of the integrals for the general case of any number of layers is not readily 
carried out unless certain approximating conditions are introduced. The field E , for the case 
of the line source and the observer in the interface (z=h=0), is conveniently written 



where 



2x Jo u +u e 
^Hl + OWOtanhwJ' (98) 



where F 2 is the wave admittance associated with the lower layers and is given by eq (48) et al. 
When hi becomes sufficiently large, it is seen that u e =U\ and the fields correspond to a homo- 
geneous half-space. 

Now, when the subsurface layers are highly conducting, such that |y™|>]>|yoL it follows 
that u m ={\ 2J ry% l )^^y m , since the important values of X are when X is of the order of |yoI. To 
this approximation 

where y e can be regarded as an effective propagation constant. It is approximately given by 

7eV (Ta/JVJ+tanh 11A 1 ( . 

T^Ll + OViVOtanh uAJ^o { } 

For a two-layer medium (h 2 -^*>), Ye=Yi/$, where 

^ = l+(7i/T 2 )tanh Tl A/ UUU; 

which is identical to eq (13). 

In analogy to eq (64) for the homogeneous half-space, the field E for the stratified medium 
can be written 

E ^ ~ l r% ' gi(7„s)<y for| 7 ^i»l. (101) 

T[y 1 -~ y 2 )x 

The factor Q 2 can be interpreted as a correction factor to account for the presence of stratifica- 
tion in the half space. A numerical discussion of Q has already been given. 

9. Magnetic Line Source Over a Stratified Medium 

In the previous section, a line source of electric current / over a stratified medium was 
discussed. The corresponding formal solution for a line source of magnetic current K is ob- 
tained by making the transformations: I->K } (T m -\-ie m o)— >iiJ, m oo, i/j, m u)-><r m -\-i6 m <ti, E Z ->H Z , 
H x =—E x , and H y -^—E y . The nature of the resulting integrals in the solution are, however, 
more complicated than in the electric-current counterpart. 

The integral form of the solution can be written, for 2<0 

H = ~ l " e ° K [ +m [e ±u ^ +n) +Bi\)e u ^- h ]e-^ x d\ (102) 

4-7T J _ m 

where the + sign is to be used for (z+h)<^0 and the — sign for (2+A)>0. The reflection 
coefficient is now given by 

W = fei' (103) 
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where K Q and Z x have their usual meaning. For a homogeneous half-space, for example, 
Zi=K x and if Aei=/i , 

it II (a; =-o — j — o — 

The coefficient R n (\) can have a pole when y\u J rylui = 0, which is the "resonance" condition 
discussed previously. This complicates the evaluation of the integral and eliminates the 
possibility of using the formulas developed for the line-current source in the previous section. 
(There was no pole in the reflection coefficient R±(\) for a homogeneous half -space, because 
there is no solution of u +Ui = 0.) 

If attention is restricted, however, to the field a large distance from the surface, or if the 
height h is large, it can be expected that contributions from the pole (which are surface waves) 
are small compared to the specular component. Then with analogy to eq (57), the field is 
approximately given by 

Ho ^Z^^ {Ko[yo[x 2 +{z+m] + R{l ^ (1Q4) 

where R u 0) is the reflection coefficient for a plane wave (parallel polarization) incident at an 
angle on the interface (2=0) where tan 0=x/(z-h). This result is only valid in the asymptotic 
sense, and hence to the same accuracy the formula reads 

H ^~^r (2) b^ +R «® wi (105) 

where only the first term of the asymptotic expansion of K has been retained. As can be 
verified by the saddle-point method (see next section), this result for H Q quoted above is the 
first term of an asymptotic expansion. 



Figure 7. Magnetic line source over a dielectric 
coated conductor. 
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10. Magnetic Line Source Over a Dielectric Coated Conductor 

As mentioned above, the integrals appearing in the solution for parallel polarization are 
complicated by the existence of poles. This occurs even in the case of a homogeneous half- 
space. Another relatively simple situation is when the upper layer is a pure dielectric and the 
region below is a metallic conductor [9]. The situation is illustrated in figure 7. 

Noting that Z 2 =0, the field for 2<0 can be expressed by 

-ie a>K r+ r ± K (^ )+ ^o-^oA eU0(z - h) l e - ixxdx (106) 



where 



Ki tanh uji\ * TJ r 

A= — - — and K x - 



Ui 



Vo (Ti + tCOCi 
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The integral has a pole where Uo+ik(fi=Q, an d tnis is denoted by \ x . Furthermore, there is a 
branch point where u =0, w T hich is when \=k . Now for small values of h x (i. e., k x Ai«l) A 
is also small (compared to one) , and therefore 



KiUih\ ifjL coh 



Vo 



Vo 



! 0-l> 



Therefore, to a first approximation, A does not depend on X the integration variable, and the 
surface impedance Z](^iy A) is almost purely inductive. The pole is then given by 



X, ^fc [l+ (ioW 2 ] 1 ' 2 ^! + (hh) 2 /2] 



(107) 



where h=h 1 (l~kl/kl) 



The integral for H is now in a form where it can be evaluated directly by the modified 
saddle-point method. A similar integral appears, in the case of a dipolar excitation. For 
the present purpose, however, it is possible to develop the first few terms of the asymptotic 
expansion in a direct manner. 




Figure 8. Integration contour in the X plane. 



The integration contour is shown in figure 8. The integration in the original integral is 
along C from — <» to +o° . The pole Xi and the branch point k can have a small but vanishing 
negative imaginary part if the conductivity of the air and the dielectric is finite. The contour 
C then passes above Xi and k but passes below — \ and — k . 

Following standard practice in contour evaluation of integrals, the contour C is closed 
by an infinite semicircle (C a +C e ) in the lower plane, with an indentation along C h , C c , and C d to 
the branch point at k . The integral along C is then equal to integral along the complete 
closed contour C minus the individual contributions of C a , C b , C C9 C d , and C e . The contributions 
from C c and C h vanish as R— >oo 9 since the real part of u ^>0. Furthermore, the integral along 
the branch lines can be obtained as an inverse series in x (which is asymptotic in nature) by 
expanding the integrand in a Taylor series and integrating term by term. The integral around 
the complete contour C is simply— 2ni times the residue at the pole \. Finally, it follows that 
the integral along C is — 2iri times the residue at \ minus the contribution along the branch 
lines. 

The final result reads, for x^>^>(h—z) 

plus terms containing -p. — r^j 77 — r^? etc. L 

r (k x) 5/2 (k x) 7/2 J 

Up= ^/Xl-k 2 Q =-ik A- (108) 



where 



The first term in the brackets is the residue at the pole X^ and the succeeding terms are from 
the branch-line integration. The latter series is asymptotic and is usable if 
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(!)<<**>. 



which requires that the pole \ x is not too near k . Such a condition is violated when the thick- 
ness hi of the dielectric coating approaches zero. It is then necessary to apply the modified 
saddle-point method [10]. 

1 1 . Fields of a Vertical Electric Dipole Over a Stratified Half-Space 

In the previous section, the excitation of the fields over a stratified half -space was by a 
line source. In this case the problem is a two-dimensional one. It is the purpose of this 
section to mention the extension to the case where the source is a dipole. This problem was 
formally solved by Sommerfeld in 1909 for the homogeneous half-space. 

A vertical electric dipole is to be located at z=—h over the stratified region of M layers. 
The dipole is to be regarded as a current element of length ds carrying a current /. 

The situation is illustrated in figure 9, where it is noted that cylindrical coordinates 
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Figure 9. Vertical electric dipole over an M-layered 
half-space. 



(p } <l>,z) are employed with the 2-axis pointing downward as usual. The primary influence in 
the upper region (z<C0) can be expressed in terms of a Hertz vector, which has only a z compo- 
nent 7rf , which hereafter is denoted \p p . Furthermore, 

Ids e~ ik ° R ° 



V= 



where 



7?o= 






(109) 



It is now assumed that the total fields can be expressed everywhere in terras of a Hertz 
vector, which also has only a z component, denoted \f/ m . (The justification for doing this is 
attained by verifying that the final solution satisfies all the boundary conditions and behaves 
properly at the source and at infinity.) For the mth layer, \[/ m satisfies 



(V 2 -7m)V'm=0, 

except at the source where 4 , a-^^ /V , for R u -^0. This suggests writing 



(110) 



(in) 
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where ^ is the "secondary influence". Since 

(V 2 __7o)^=0, (n2) 

then 

(V 2 -7o)^ s =0. (113) 

In cylindrical coordinates and for azimuthal symmetry, the equation for ty m becomes 

Gli'li+S-^-)*-" ' (114) 

It can be readily seen that solutions of this equation are made up of linear combinations of 

* w ™*Jo(Ap) 

r (xp) 

where, as usual, u m =(\ 2 +yl) H , and 7ro= (^m/%«— € w /x m co 2 ), for m = 0, 1, 2 . . ., Af. Again by 
convention, the real part of u m is to be taken positively. Since the resulting solution is to be 
finite on the axis p=0 when z^—h, the F (Xp) Bessel function can be rejected because it 
behaves as log p for small values of p. The J (Xp) Bessel function on the other hand, is finite 
at p=0. 

As in the previous solution for the line source, the parameter X can take any value (with 
a few exceptions) . The general solution for \p m is then written in the form 

^m=J o a5 K(X)6-^ 2 +6 M (X)^- 2 ]Jo(Xp)dX, (115) 

for ra=l, 2 . . ., M— 1. In region (0), for 2>0 

*>=**+ f &o(X)6 W0 *J (Xp)dX, (116) 

and in region (M), for z^>z M -\ 

lfc*=f a M (S)e' u ^dk. (117) 

To express \p p as an integral of the appropriate type, use is made of the Fourier-Bessel 
transforms, which read 



and 



f(p)=£ gWJoM*d\, (118) 

^(X)=J o <X> /(p)Jo(Xp)prf P , (119) 



subject to the conditions that the integrals exist. Letting f(p) = e~ ik ° p /p, 

g(X)=J «-***> J (\p)dp. (120) 

This integral is a standard type [5] and is given by 

«W"o?zJBJrv (12D 

and therefore 

e~ ik op _ r°° J (Xp) ^^_e~ ik ° R ° 
P Jo 



\d\=- 



p Jo (\ 2 -H) H Ro 
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(122) 



To generalize the preceding to any other value of z, the integrand must contain the factor 
e *u {z+h) b ecause e~ ik()R() /Ris a solution of the wave equation and necessarily becomes e~ ik ° p /p 
as z-^—h. Consequent^, 

Jn(\o)f ±UQ{Z+h) 

-\d\, (123) 



Bo Jo u 

where the sign of the exponent is chosen so that the integral converges, 
(i. e., + for {z+h)<0 and - for (z+h)>0). 

The z component of the Hertz vector for the upper region (z<CJS) can then be written 



HT[ a 



«o L - +6(X)— — U(Xp)Xd\, (124) 

where 

a = -r- -. — and b(\) 
47rtcoe 

is to be determined from the boundary conditions. Now 

J&¥». jp _n w _/_ -.2 , ^ 

and 



^ mp== d~^z } ^™* =0 > E mz =( — 7w+j^p )^ ( 125 ) 



ff TOP =0, fl^-jjA^a, fl^=0, (126) 

^/x w co Op 

and since 2? OTp and fiT*^ are continuous at the interfaces z m (m=0 to M—l), it follows that 
tymfbz and (o- OT +iwe TO )^ m are also continuous at the same interfaces. This leads to 2(M— 1) 
equations to solve for the 2(M— 1) unknown coefficients. It readily follows that 

— =«(X)=^ (X)+Zi(x) > (127) 

where 

t- />\ 7^ a\ Z W+1 (X)+7T W (X) tanh^/i^ noQ \ 

Z m (X)-K,„(X) • Km(}i)+Zm+l(K)t&nhUnhm - (1-8) 

The algebraic form of 11 (X) is identical to eq (5) for the plane-wave reflection coefficient for 
parallel incidence. It is thus clear that the field of a dipole over a stratified half-space can be 
regarded as a spectrum of plane waves whose angle of incidence and reflection is related to the 
variable X by \=&o sin 6. In this case, the wave normals generate a family of cones coaxial 
with the 2-axis. 

Equations (127) and (128) in conjunction with eq (124) constitute the complete formal 
solution. Approximate evaluations [10] can be carried out by using the same approach as in 
the previous sections or by a modified saddle-point method [11]. In view of the similarity 
with the line-source excitation, it does not seem justified to discuss this in detail, and the 
essential results will simply be quoted. 

As has been demonstrated, Z x {\) for a highly conducting half-space or even for a dielectric 
coated conductor, is mainly determined by the electrical properties of the layers. In other 
words, 1\{\) can be replaced by 

^-^ g) 

v u +ik A 

where A=Zi/r) =Z 1 ]120ir is assumed to be constant. The form of ^ is then 

+o=t^- P [ e ±u iz + n )+ ^-lhA e u iz-H)lj MdKm (130 ) 

4irlueoJo L Uo + ifaA J 
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This can be rewritten as follows: 






where 
and 



P -Jo (« +i*oA)«o J ° (Xp)rfX ' (132) 



This integral has been previously evaluated [10, 11]. Subject to |A|<<1, it is given by 

e -ik Ri 

where 

w=p ( 



P^i(wp 6 )^ e~ w ertc(iw*) —5— (133) 



and 

erfc (ot j ^) ==— 1/o e~ x2 dx- 

-y7T Jiw /2 

The expression for P can also be written 

P=f [l-W]S^> (134) 

where 

F(w) = l-i(Tw)V w erfc (w 4 ). (135) 



It is noted that when z=h=0, 



. Ids e~ ik cp v . 



and the vertical electric field is given by 

E 0l J^-e-^F( Pe ) ioxk oP »l. 

The function F(p e ) can be regarded as the correction to the field of a dipole on the surface 
of a perfectly conducting plane. For |^|<C<^1, it approaches unity. F(p e ) has the same func- 
tional form as the ground-wave attenuation of Norton [6] who presented numerical values for the 
case where (in the present notation) b is in the quadrant 0° to —180°. For a homogeneous 
ground A has a phase angle in the range 0° (for a perfect dielectric) to 45° (for a good conductor) . 
The corresponding values of b are —90° and 0°, respectively. For a stratified ground, however, 
the phase angle of A may be outside this range. In the particular case of a two-layer conducting 
half-space, the parameter p e can be written 

Pe=PiQ\ where p x = — ^ (V^o) 2 ' 
and where 

rt-lrtl,«t~ ( 7i/7 2 ) + tanh7iAi 

V-\V\ e = l + ( 7l /7 2 ) tanhTxAi' 

in terms of the propagation constants 7l and y 2 of the upper and lower layers and the thickness 
hi of the upper layer. This function Q is discussed in section 2. From figure 3b it is seen that 
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FIGURE 10a. Amplitude of propagation factor for an 
impedance boundary. 
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Ficure 10b. Phase of the propagation factor for an 
impedance boundary. 



for a highly conducting substratum q can become positive. Consequently, the phase angle 6 
can exceed and would approach 90° for IY2/Y1I-* 00 , and TiAi<<1. 

In analogy to the terminology introduced by Sommerfeld for the homogeneous half-space, 
p e is here described as an effective numerical distance. The propagation factor Fe^lFele'*®, 
is shown plotted in figures 10a and 10b, as a function of \p e \ for a range of b values. The curves 
for 6<0 correspond to those computed by Norton, and it is interesting to note that for this 
region, \F e \ never exceeds unity. On the other hand, when 6>0, \F e \ may exceed unity, and 
for 6 = 90° this effect is most pronounced. Apparently, in the case where b is positive, the 
energy is being guided to some extent along the surface. This effect can be seen in the asymp- 
totic development of F(p e ), which, for be|>>l, reads 



F(Pe)^- 



1 1-3 1-3-5 1-3-5-7 



when — 27r<6<0, and 



2p e (2p e )> (2p e )* (2p e y 

1 1-3 1-3-5 



F(p e )^-2i^7rp e e- p 



when 2tt>J>0- 



2p e (2p e y (2p e y 



(136) 
(137) 



The term —2i^wp e e~ Pe has all the characteristics of a surface wave. It is not present in th e 
asymptotic development when b is negative. At 6 = 0, 27r, or — 2t, p e is real, and this term 
vanishes asymptotically. The trapped surface wave is most predominant when 6 = 90°, 
which corresponds to A or Z x being positive imaginary. For example, if the surface is coated by 
a thin dielectric film of thickness hi with a dielectric constant ei, 



Z x ^ i^h x [ 1 — €b/ei] ^ iX l , 
which is purely imaginary. Then, since 

— ik p 



Pe 



(Zi/W 2 , 



(138) 



(139) 



Pe=\p e \e iT/2 , or b = ir/2 radians=90°. The attenuation factor now has the asymptotic 
development 



F(p e )^-2e iST/A ^Tr\p e \e- i 



\Pe\- 



3 



15i 



2\Pe\ ' ±\Pe\ 2 »\Pe\ 



(140) 
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where 

"2 r,l 
and the vertical electric field is 



Pe\- 



E*eJ^.e- a +Ffa), (141) 



which shows that the (trapped) surface wave component varies as 

k p 



which has the characteristics of a cylindrical wave traveling in the positive p with a phase 
velocity 



o+sr 



times that of light. 

The (trapped) surface wave is not excited when b is negative, and then the distant electric 
field varies asymptotically as 

(M 2 

which has a phase velocity equal to that of light. 

When the source dipole and the observer are both raised above the surface (i.e., h^>0 and 
— s>0), the field in the upper half-space can be expressed conveniently in the asymptotic 
sense as the sum of three partial fields in the manner 



where e=0 for 6<0, and €=1 for 6>0. An asymptotic development shows that 

,-ikR Q ^ e ~ikj 

~R^ + C+A~R l 



e -iJcR C-A e~ ikRi ,. Ah . 



{1 1-3 I.3.5 ^ e -ikRi 

^(l + C7A) 3+ 2p 2 (l + f7/A) 3+ 4^ 3 (l + C/A) 7+ • * 7 ~Wf (144) 

and 

^-^,[2^7^-] -^-, (145) 

where 

C=— 5— = — ^-J, fc=:27r/wave]ength=fc , 

Xii ill 

and A=Zi/riQ. 

In terms of ^ , the vertical electric field components for & p>>l is given by E^^klil — 

In the above, \p a can be identified as a geometrical optical term being just the primary 
field 6~ ifci2 o/R and a specularly reflected component e~ i1cR i/Ri modified by a reflection coefficient 
(C—A)/(C+A). yp b is an asymptotic expansion containing terms varying as R{~ 2 , Rr 3 , etc. 
\f/ 3 , which is only present for positive b values, represents a (trapped) surface wave and has a 
phase velocity less than that of light. When b is negative there is no (trapped) surface wave 
present. 

Sometimes \p a itself is called a space wave and \j/ b a surface wave. This usage would cor- 
respond to that of Norton [3]. This would be an obvious designation when the half -space is 
homogeneous because \f/ $ is zero, and there is no trapped surface wave excited. 

Equations (143), (144), and (145) enable one to compute the field in the upper half- 
space (?<0) for an electric dipole located at z=—h in terms of the ratio A of the surface im- 
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pedance Z x to that of free space 770. As a first approximation, Z x could be replaced by Zi(X) 
for \=0, which would be the normal surface impedance. Clearly, a better value would be 
Zi(\ s ), where X s is the saddle point of the integral in eq (130) or (132). This would mean that 



A= j 

Vo 



where \ s =k S, where S=p/R 1 =- y /l- 



-C 2 . For a homogeneous half-space, with Mi=A<o, 



Vo 






2 \H 



(146) 

(147) 
(148) 



The coefficient (C—A)/(C+A) in eq (143) now can be identified as the Fresnel reflection co- 
efficient for a plane wave with parallel polarization incident at angle 6. The equations are 
now in complete agreement with the results of Norton [6] for a vertical electric dipole over a 
homogeneous conducting ground. 

The corresponding treatment for a vertical magnetic dipole located at z=—h over the M 
layered half-space is almost identical to the above. The quantity ^ is then to be identified 
with the z component of the magnetic Hertz vector, and then A is to be defined by A=r7 Fi, 
where Y Y is the surface admittance at the interface 2=0. 

12. Appendix 
12.1. Effect of a Transition Region on Reflection 

The reflection problems considered above were entirely concerned with layered media. 
Each layer was assumed to be homogeneous and bounded by a sharp interface. 

In certain instances the dielectric constant or magnetic permeability can be regarded to 
vary in a simple analytical fashion. In such cases, it is sometimes possible to obtain a solu- 
tion of Maxwell's equations in terms of known functions. The reflection coefficient for a 
medium in which the dielectric constant may be expressed as a power of one coordinate has 
been considered in great detail by Forsterling [12]. The hyperbolic transition has been con- 
sidered by Rydbeck [14]. 

In this appendix, a treatment of the linear transition between two lossless dielectrics will 
be given. Actually, this is a special case of Forsterling's model, but the mathematics is some- 
what simpler. 

The situation is illustrated in figure 11. A plane wave is incident at angle from the air 
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Figure 11. Illustrating a plane-wave incident on a 
linear transition region. 
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to an interface at 2=0. The electric vector is perpendicular to the plane of incidence. The 
dielectric constant between z=0 and z=d is a function of the coordinate z and denoted €i(z). 
The region below, z^d, has a dielectric constant e 2 , which does not vary with z. The magnetic 
permeability of the whole space is taken as /z and assumed to be constant. 
The electric field, which has only a y component, can then be written 

2£ = [ e ~ifoz cos 0_|_ ]$ e ihz cos 0ol e ~ikox sin do f or £<^(), (149 ) 

and 

E 2= T e -ik2 2 c:se2 e -ikoxsme a f or 2>d, (150) 

where 6 2 is related to di by k 2 sin 2 = =&o sin O - The constants .R and T are to be determined by 
the boundary conditions at z=0 and z=d. 

Within the transition layer, the field components satisfy Maxwell's equations 

~^-==iwuH lxJ (151a) 

iix o)H lz =ik sin d E u (151b) 

and 

=^=i[t, (2) - 6 sin 2 O ]«E„ (151c) 

oz 

where E { is the y component of the electric field in the transition region. For a linear transi- 
tion. 

ei (2) =«o(l +02), (152) 

where /3 = (« 2 — « ) /e d. A new variahle IF is now introduced, defined by 

TJ7 TI7 /1 I # 2 V /2 1 TI/ 2 (''l C0S3 •) /I M\ 

W= ^( 1+ co^) and W » = 30 (153) 

Equation (151a) and (151b) can then be written 

<«§=«, and y^^ia&Ev (154) 

where e=e l (z) — e sin 2 0. (155) 

These equations are combined to give 

(^+imw+^) El=0 - (156) 

But, since e= const. XW m , 



(157) 



_1_ (U_ J_ 
2edW~SW 

Furthermore, if the substitution 

E X =W^ (158) 

is made, ty satisfies the equation 



w§-^+(^'-hy°- 



(159) 



which is Bessel's equation of order 1/3. Independent solutions are Jj^(«ju' W) and J.^ 
(w/x M W). 
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The general solution for E 1 is then written 

£\=(I^/TF ) H [A/ H (c^ sin O ), (160) 

where A and B are constants. The magnetic-field components are obtained from 

(161) 



IIi x =-. 5—= and H u =— —E u 



ijJL CO dz 



Mow 



and need not be written out. 

The four unknown constants R, T, A, and B are obtained by matching the tangential 
fields at 2=0 and z=d. That is, 



E =E, ' 

dEo = dEi 
dz dz 



E X =E 2 " 

and m^m 

-. {) dz dz 



1 1 r,i>) 



Omitting algebraic details, the expression for the reflection coefficient can be finally written 

1-6? 



R= 



: 1+G' 



where 



vhon 



with 



and 



(boq-i+aopj—iiaobx—aibo) ' 



6 =«/-H( a; o), qo=J-n(xo), 
th^Jiitedf 2>i=Jy 3 (xi), 
b x = J^ A {x l ), q l = J. % (x ] ), 



(163) 



(164) 



0*0= 



"iV 2 -l 



4tt aN* 4?r 



'N 2 -l 3 



COS o =;r--eOS 0o> 



tf= 



(e 2 /e o -sin 2 0o)" 2 



cos0 

Employing tables of fractional order Bessel functions [15], the magnitude of the reflection 
coefficient \R\, is plotted in figure 12 as a function of a or (d/\ )cos O . The abscissa is propor- 
tional to the thickness of the transition region. Two values of N are chosen, 3 and 5. 

It is immediately apparent that the effect of the transition region is to reduce the reflection 
coefficient. For small values of (d cos O /X O ), the reflection is governed almost entirely by the 
(normalized) relative refractive index N. In fact, as d— > 0, it is readily seen that 



#= 



1— N_cos d — (e 2 /ei — sin 2 o ) K> 
''1+N~ eos0+(e 2 /ei-sin 2 0o) 1/2 ' 



(165) 



which is the classical Fresnel formula for perpendicular incidence. As seen in figure 12 \R\ = 
0.5 for 2V=3, and d=0; \R\ = 0.666 for N=5; and d=0. 

In the language of transmission line theory, the transition layer is acting as a matching 
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Figure 12. Reflection coefficient for a linear transition layer. 



section. Its function is to match the impedance in the air region (z<C.O) to that of the lower 
region (z^>d). 

The simplicity of the above analysis is due primarily to the choice of the polarization of 
the incident wave. In the corresponding problem for parallel incidence, the resulting equations 
for the magnetic field component in the transition la3^er are not expressible in terms of Bessel 
functions. 
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